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' The continuum Faddeev equations for the neutron-neutron-alpha (n+n+a) system are formu- 

, lated for a general interaction as well as for finite rank forces. In addition, the capture process 

■ n+n+a — >^He+7 is derived. 

> ! 

Q PACS numbers: 21.45.-v, 25.10.+S, 25.40.Lw 

; I. INTRODUCTION 

(N : 

I— I In recent years the study of quantum halo systems experienced increased interest in the nuclear as well as the atomic 

■ few-body community. For a recent review see Ref. The nucleus ^He is of particular interest, since it constitutes 
' , the lightest two-neutron halo nucleus with a ''He core. Being an effective three-body system, the properties of the 

Q ' ground state have been explored using either Faddeev 0-01 or Hyper-spherical Harmonics (HH) . More recently 
^ ] the ground state has also been calculated with multi-cluster methods [9l-ll7t as well as with GFMC ^Sj] . Those multi- 
, ^ , ■ cluster methods include various techniques like the microscopic dynamical multi-configuration three-cluster model [§] , 
the stochastic variational method [Tol|, the multi-cluster dynamic model (MDMP and AMDMP), the hybrid- TV model, 
psj a combination of the cluster orbital shell model (COSM) 15] and the extended cluster model (ECM) JJJ, the refined 
^ resonatin g gr oup method [ill , [l^ l , and the coupled-rearrangement-channel variational method with Gaussian basis 
functions [175. In. addition, the beta decay to the alpha+d continuum has been studied [l^, [l^l. Of interest is also the 
two-neutron capture process '*He(2n,7)^He as a possible route bridging the instability gap at nuclear mass A=5 pTj . 

The situation is quite different in the continuum of two neutrons and an a -p article. There is a well established 2"*" 
resonance [l^], but further resonant structures are still under debate (tI. [23l - l30| . Up to now more indirect approaches 
in understanding the resonance structure have been carried out, e.g. the four-body distorted wave approach, leading 
to three-body continuum excitations of two- neutron Borromcan halo nuclei (23l. [28j. Furthermore, complex scaling 
in Coulomb break-up reactions has been employed [31.] . In addition, an extension of the HH method on a Lagrange 
mesh [30| has been used to study three-body continuum states. This is at least a four-body problem with great 
uncertainties about the reaction mechanisms and the interactions entering these much more complicated systems. 

Thus, the currently predominant approach to continuum calculations for the pure n + n + a system is the HH 
method [tI. IsI. [26l - l28l . l32j] . A Faddeev approach is to the best of our knowledge still missing. Only for the ^Li nucleus, 
] a Faddeev treatment of the deuteron-alpha (d-a) system has been employed [33j , which howeverjdid not have to face 
the challenge of three-to-three scattering. This also refers to the pioneering work by Koike [13, [l^ on d-a scattering. 

The aim of this investigation is to fill that gap. For the n + n + a system one faces the situation of three free particles 
being in the initial channel and leading again to three free particles in the exit channel. In other words, one has to 
deal with three-to-three scatte ring . Scattering of three free incoming particles to three free outgoing ones in a Faddeev 
approach has been initiated in [Sa 113] in the context of the three-body photo disintegration of ^He. This path has also 
been followed in the same context and in a Faddeev approach by Meijgaard and Tjon [381] . Into the matrix element for 
the photo disintegration enters the three-nucleon to three-nucleon scattering wave function, which has been evaluated 
in Ref. [38] and then inserted into the photodisintegration matrix element. However, evaluating the wave function 
is a completely unnecessary complication, since this process is initiated by the three-nucleon bound state. One can 
directly derive a Faddeev equation for the three-body break up amplitude, in which the driving term contains the 
action of the current operator on the ■^He ground state. Then, the complete final state interaction is generated by 
a Faddeev integral kernel for the amplitude. This considerably simplifies the technical part of a calculation, since 
no disconnected processes occur. This very procedure has been pioneered in Refs. [36l \3l\ and is being applied in 
state-of-the-art calculations, see e.g. [sollioj. 

The same procedure can trivially be adapted to the capture process n + n + a ^•'^He+7, as will be displayed in the 
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present investigation. This capture process is relevant for the production rate of ^He in astrophysical environments 41 1 
characterized by high neutron and alpha densities e.g. those related to supernova shock fronts. In e.g. Ref. [2l| this 
three-body process is approximated by sequential two-body processes, whereas in principle a genuine three-body 
reactions needs to be calculated. Furthermore, the nna — > nna amplitude is relevant for determining the the next 
order coefficient [i^l in the virial equation of state in low-density matter [4^ . 

From a technical point of view the Faddeev approach to the n + n + a continuum is strongly needed, since the 
currently predominant approach, namely the Hyper-spherical Harmonics (HH) approach, still faces open challenges. It 
is already known that in the break up process n + d ^ n + n+p a, strong FSI peak appears for the n — n subsystem. In 
the Faddeev approach using Jacobi momenta this can be mapped out correctly, whereas when changing to the hyper- 
spherical angle, the convergences is quite poor for this particular configuration. In the HH method, the control of the 
coupling potentials can be a painful exercise, whereas in the Faddeev approach using Jacobi variables the dynamics 
is perfectly well under control in all details. This same situation must be expected in n + n + a scattering, where 
the three-body S-matrix is characterized by continuous quantum numbers describing how the energy is distributed 
over the relative motion. There are strong initial and final state interaction peaks, which in a discrete representation 
through hyper-spherical K-quantum numbers are difficult to map out correctly. As stated above only a technically 
reliable approach as the Faddeev one will guarantee the validity of the results when searching for ^He resonances. 

The paper is organized as follows. In Section II we derive the coupled Faddeev equations for the three-to-three 
scattering amplitude, followed by a partial wave decomposition in Section HI. The Faddeev equations will be solved 
by iteration yielding a multiple scattering series. This will be outlined in Section IV. Since most of the Faddeev type 
investigations of the n + n + a system are based on finite rank forces, we also present in Section V a continuum 
formulation based on separable forces. Furthermore, we discuss the the unitarity relation for the three-to-three 
amplitude in Section VI. Finally the capture process n + n + a ^•^He+7 will be discussed for the Faddeev scheme 
in Section VII. Then we summarize in section VII. Technical details about the partial wave decomposition and an 
efficient way of treating the three-body singularities are given in the Appendices. 



II. THE FADDEEV EQUATIONS FOR THE NNa SYSTEM 



In developing the formal expression for the transition amplitude between three free particles interacting with 
short-range, strong interactions, we start from the triad of Lippmann-Schwinger (LS) equations (4^ . [45l | acting on a 
three-particle initial state given by 

$i+^ = |pJ^+^|q«>, (2.1) 

where Ipq)'^-* is a two-body scattering state, and the index a = 1, 2, 3 indicates the three choices of pairs characterized 
by the third particle, the spectator. Furthermore, — J^p^a ' where Vg {(3 — 1, 2, 3) are the pair forces. Three- 
body forces can in principle be incorporated in a straightforward fashion. However, we will only concentrate on 
two-body forces here. The triad of LS equations, 

+ G„ (2.2) 

define the scattering wave uniquely. The channel Green's function is given by G"^ ^ {E + ie — Hq — Va)^^. We use 
standard Jacobi momenta Pa and and their quantum numbers as basis states. 

By suitable multiplication of the three equations in the triad from the left by Vp one obtains the transition operators 

Uao = {Vp + V^)^^+\ with 13^ ^ a, which fulfill the set of equations 

Uc.0 = + J2 t/pGoUpo, (2.3) 

where $o = |p)|q) is the free three-particle state. 
The three-body break up operator is given by 

C/oo^5]l^^*^+^ (2.4) 

7 

Again, from the triad follows 
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Iterating Eq. (|2.3p one obtains the multiple scattering series 

Uoo = *^*o + ^7^0 ^/^^o + Y ^7^0 Y Y ^"*o + ■ ■ • • (2.6) 

Instead of working with the coupled set of Eq. (|2.3p and the relation of Eq. (|2.5I) for the three-body break-up 
operator, one can generate the multiple scattering series directly by decomposing Uqq as 

[/oo = 5]f/7, (2.7) 

7 

and choosing to obey the coupled set of Faddeev equations 

Uj = t^+ t^Go Y Ua- (2.8) 
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Indeed, iterating Eq. (|2.8p and inserting the result into Eq. ()2.7p leads exactly to the multiple scattering series from 
above. Explicitly, we have a set of three coupled equations 

Ui = h+hGo{U2+U3) 

Us = h + t3Go{Ui + U2). (2.9) 
We also observe, that comparing Eqs. (|2.7p and (|2.4p leads to 



EE (2.10) 

When considering the n-|-n-|-a system, we need to incorporate the identity of the two neutrons. Fixing arbitrarily 
the alpha particle as spectator and label it as "1" and the two neutrons as particles "2" and "3" , the scattering wave 
function ^q^' must be antisymmetric under the exchange of particles "2" and "3" . Thus, defining the transposition 
operator P23, the scattering wave function must fulfill ^23^0^^ — ~^'o^^- Using this in Eq. (|2.10p leads to 

U3 = -P23C/2. (2.11) 

Thus, for the n-|-n-|-a system we only have 2 coupled equations, 

Ui = h+hGo (I-F23) U2 

U2 = h+ i2Go hP23U2 + Ui) (2.12) 

More precisely, one has to apply the driving terms to the free state $o,a, which is antisymmetric under exchange of 
the two neutrons: 

$0,a= (l-P23)|Piqi)|0TO2TO3)|0~), (2.13) 



leading to 



Ul^O.a = tl$0.a+ilGo(l-P23)C/2«'o.a 

f/2$0,a = t2%,a+t2Go{-P23U2%,a + Ui<i>o.a)- (2.14) 



The full break-up operator is given by 

i7oo$0,a = f/l*0,a + (1 " P23)U2<S?0,a- (2.15) 

For the on shell break-up amplitude one has to evaluate the matrix element ($q a\Uoo,a\^o,a), where in the final state 
momenta as well as spin magnetic quantum numbers are changed, 

'i>o,a = (l-^23)|p'iq'i)|0m>^)|oii). (2.16) 
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III. PARTIAL WAVE DECOMPOSITION 

In order to solve the coupled equations, Eqs. (|2.14p . two sets of partial wave basis states are needed: 

1 1^ 



22 



\piqiai) = C(jiAi J, ^lAfi - /ii)|pi(?iSi)ji/ii)|giAiMi - /ii 

Ml 

\p2q2Ct2) = '^C(j2hJ,lJ'2M2 - /X2)|p2(^2S2)j2Ai2)k2(A2-)/2M2 " M2) 



M2 

'1 1 



The details of a partial wave decomposition of Eqs. (|2.14p is well known (see e.g. [1^), and we refer to j40| for details. 
Employing the states of Eq. (|3.1I) . the coupled equations, Eq. (I2.14p read 

^^^^^ Mb'iPi,i?,J 

+ + j dx j dq'^q'^^ t^,^{p[7r,{q[q'2x),Eg,) Go{7:i{q[q'2x),q[) 

^Ga'^a'^{q[q2x){TT2{q[q2x)q'2a'^\U2,a 

{pW2<^2\U2,a = 

12 

92 

+ dx dq[qj^ ta'^{p'2TT3{q2q'ix),Eq'J Go{Tr3{qWix),q'2) 



^Ha'^a[iqWi^)Mq'2q'ix)<iWi\Ui.a 



where 



with Ai = 2Ai + 1. 



dx J dq^ q2 ^ ia^(P2'^5('?2'?2 x),Eq,J Go(7r5(g^92 2^)><?2) 
Xl^a^4"('?292"2;)(7r6(g292"a^)92""2V2,a , (3.2) 



CZ'+""'{0i) = (l + (-)'^+^0 (^^Si,m2m3W(ZiSiji,mi,,m2+m3) 

^ ^ mil 

(jiAiJ,mi, +m2+m3,0,M)y,,™ (0i,O)W^ , (3.3) 

V 47r 



D^-™^(0i)^i?™-'"^(0p,0,J 

.1/ 



'^siiXl(-^'2^2>^''^'^^'"^) M2 2J2'^-"^3,"l3j y,^^_„3(p2) 

>''2t:I2,M' ~ 1112,1712) Yy ,,i,_^_.^^{q2) (3.4) 



and 



^41 E U2I2J',^^M' - fi) ( l'2h'2,fi- m2,m2 ) r4^_„,(4) 



2-' 
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The 'shifted' momenta tt^ are given as 



-^2^^2:^^'-A*-"l3,m3 ) ^X^M'-p-mafe) • (3-5) 



71"! 

7r4 





-1— r\n^ /t' t* 




+ 2/?gi92a^ 










^^2^ + 92" 






2 + 2^g^q2"'x 



= \/92 +/3 92 ' + 2/3g^'Z2"'x, (3.6) 



where 



a 
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1m + rriQ, 

2(to + Tfla) 



-P - ^ ■ (3.7) 

TO + TOq 

Here to is the neutron mass, and TOq the mass of the "^He nucleus. 

We refer to the Appendix |X] for some details of the derivation and the expressions of the purely geometric quan- 
tities G a'^a'^{q'iq2x) , -ffa^a'j (9291^)' ^^"^ ^a' a" ^^2l2 ^) ■ Furthermore, ta'^(p[pi, Eq-^) is the two-neutron f-matrix and 
ta'^{p2P2, Eq^) the One for the neutron-a pair. 

Due to the free Green's functions Gq and the x- integration over it, one encounters the well known logarithmic 
singularities of any three-body problem. These singularities can be reliably treated [1^ H^l- However, the method 
suggested in Refs. [3 11^ j appears to be beneficial for here, since not only kernels contain logarithmic singularities 
but also the driving terms. We illustrate this new method with an example in Appendix [Bl 

IV. THE MULTIPLE SCATTERING SERIES 

A well established way to solve a coupled set of Faddeev equations is to generate the multiple scattering series. For 
the 3N system is is laid out e.g. in Ref. [s^. Schematically the Eqs. p.2p have the form 



U = U^°'>+KU (4.1) 

which, when iterated yield 

U = C/(") + iJ^i) + C/(2) + ••• , (4.2) 

with 

C/(") ^KU^"~^\n^ 1,2,- ■■ . (4.3) 

The first few terms of this series are depicted in Fig. [T]The driving terms of Eqs. p.2p . sketched in the upper row of 
Fig. [Hare necessarily disconnected, since a two-body t-matrix can not act on three particles. 

Let us consider the terms of second order in the two-body t-matrix (indicated in the second row of Fig. [IJ: 

{p[q[a'M'l^ 

(l + (-l)''i+^'i) fdx f dq'^qi t^,^{p'^TT^{q[q'2x),Eq,^) Go{TTi{q[q'2x),q[) 



12 ^ 
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Ga[a'^{q[q2x) ta'^{Tr2{q[q2x)p2,Eg^) Z?"^" (6*1 ) 
ta[{p'i'^l{q'iq2x),Eq'^) GQ{TTi{q[q2x),q[) Ga[a'^{q[q2x) 

t^'^{7T2{q[q2x)p2,E^,)D'^,^^"'%ei) 



(4.4) 



The only singular function under the x-integral is the free Green's function, which leads in the q'l — q2 and q'l — q2 
planes of external momenta to the well known logarithmic singularities. The same is true for 

(P292«2|C/i^ ^ 

dx / dq[q^^ta'AP2T^3iq2qix),Eg' ) Go{7r3{q2q[x),q'2) 



Y^H^'^a'MWlx) ^i^L^ t^,{n4{q'2q'ix)pi,Eg,) 

- f dx f dq2 ({2'^ ta'^{p'^T:Aq2q2 x),Eq,_^) Gfi{T:Aq'^q2 x),q'^^ 
^iq'2 - 92) 



6{q2 - 92) 



92 



^.'"^(9^2 x)p2,Eg,)DyrA9^) 



^2 — t.'''iMq2q2 x)P2,Eg,) D^r^e,) 
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dx ta'^{p'2Tr3{q2qix),Eq'J Go (713(172 91 •■z;), 92) 

(92912;) t^-iMqMpi^E,,) c:^^^+"'Ae,) 



— dx 



ta'^ ip'2'^5 (9292a;) , Eq^ ) Go (tTs (52922;) , (?2 ) 

E^a;<' (92922;) t .,(7r6(g^g22;)p2,i?,J oyr^e^) 
"2" 

ta'^{P2T^5{q2q2x), Eq;^) Gq (tTj ((72922;) , ^2 ) 

E^a;4" (92922;) t .,(7r6(g^g22;)p2,i?,^J ^;;-™^(0i) 



(4.5) 



Indeed, the free Green's functions Go lead in the only remaining x-integral to logarithmic singularities. 

In order to safely apply the kernel to the previous amplitude, that amplitude has to be a smooth function. This is 
only the case for the next higher order, being of third order in t, sketched in the third row of Fig. [TJ The third order 
term reads 

{p[q[a[\u['^l EE 

'l + i-lf^+^'A j dx j dq'2q'2^t^>{p[7ri{qW2x),Eq,)Go{MqW2x),q[) 



E Ga'^a'^iq[q2x) 



dy ta'{'^2{qW2^)'^3{q2qiy),Eq') Go(7r3(g29l?/),92) 



E^«i«;'(929iy) t„;'(^4(929i2/)pi,i?,J G^?+"^(0i) 



dy 



ta' (71-2 (9l 922;)71'5 (9292?/) , Eq> ) Go (tTs (9292?/) , 92) 
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' " 2 2 2 2 

"2" 
"2" 

Correspondingly one obtains 

{p'2qW2\u^2l ^ 

^ j dx j dq[q[^ ta,^{p'^7r3{q!2q[x),Eq:J Go{n3{q!2q[x),q2) 

y iai(7r4(q2'7'ia;)7ri(gi(722/),£',j) Go(7ri(qiq2y),g'i) 

"2 

- ta'^{T:A{q'2q'ix)'Ki{q[q2y),Eq:^) GQ{'Ki{q[q2y),q'i) 

Y.Gc.',c.'S<l'li2y) to.'S^2{q[q2y)p2,E^,) by^'^^'^ie^) 



(4.6) 



dx I dq^ {q2 f ta' (P27r5('?292 x),Eq- ) Go(7r5(q2'72 2;),(?2) 



dy t^''' {neiq'^q^ x)7r3{q2 qiy),£',^) Go(7r3((72 '7iy),92 ) 



E^a'a"' (9292" 



m2+m3 / 



E-^4"«i('?2"<7iy) to,'^{'^i{q2 qiy)pi,Eq^) G™^ 



^a"' {T^Q{^2q2 x)'r^5{q2 q2y),E^,,,) Go{TT5{q'2'q2y),q'2') 



dy 



' (92 92^) (7r6(g2 '72y)P2, £^,2) ^" 



1) 



^q"' (^6(9292 a;)7r5(g2 q2y),E^n,) Goinsiq^ q2y),q2 ) 



E (92"92y) {Mq2q2y)P2, Eg,) 



(4.7) 



All three-fold integrals in Eqs. (|4.6p and (14.71) are of the same type: two angular integrations, where each one leads 
to logarithmic singularities in the corresponding momenta, one of which is external and the other the intermediate 
integration variable. It is not difficult to see that the intermediate momentum integration over products of logarithms 
leads to smooth functions in the external momenta. Therefore the third order amplitudes in t can serve as driving 
terms for the application of the kernels, and thus leading to all higher order amplitudes: 

{p'.qWM^l = 

(^1 + (-l)'i+''i^ j dx j dq'^q^^ ta,ip[TTi{q[q'2x),Eq'J GoiTTiiq'iq'2x),qi) 

^Ga'^a',{q[q'2x) {TT2iqW2x)q'2a'2\U2"a'^^ 

{pW2a'2\ut^ = 

J dx J dq[q[^ ta,',{p'^Tr3{q'^q[x),Eg,J Go{Tr3{q2q'ix),q'2) 
J2h.wMWix) {nMq'ix)q[a[\u[l-'^ 
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XI 4" ('^2 {'^6{q2Q'2 X)q'2 a'2 \U2'^a^^^ , (4.8) 

/// 

with n = 3, 4, • • • . 

The resulting series J2'^=3(pWi'^'iWi^a ^"^d J2^=3(pW2'^'2W2"a can safely be summed e.g. via Pade summation. 
For the corresponding three-nucleon amplitudes the above considerations were made in Ref. (ssj . 



V. FINITE RANK FORCES 

So far, Faddeev type studies of light nuclei treating the discrete structures have been based on finite rank forces [33|. 
Therefore, it appears useful to also formulate the nna system in the continuum in this fashion. For the sake of a 
simple notation we choose a rank-1 separable t- matrix, 

ta{PP',Eq) = ha{p)Ta{q)ha{p') ■ (5.1) 

Then Eqs. ((3?2]) take the form 

{pWia[\Ui^a = 

+ (l + (-l)''i+^9 y dxyd<7iq;2^V,(yi)Ta; (£;,;) V,(7ri(q;<Z2a^)) 

Go{TTl{q[q2x),q[) G a'^a'^{q'iq'2x) (7r2(gi(J22;)g2"2|C^2,a 

= v,(p'i)^oi(?;), (5.2) 

where the new unknown single variable amplitude is 

Za.'Sq'i) = 

^M_^ (£;,JV,(pi) C™;+™^(0i) 

+ (l + (-l)''i+^'i) JdxJ dq'2q2'Y.T^'^{E,,)h^,{MqW2x)) 

"2 

Go{TTl{q'iq2x),q'l) Ga[ai,{qW2x) (7r2 (gig22;)g2"2 I C^2,a . (5.3) 

Similarly, the second equation, Eq. (|3.2p becomes 

(P292"2|t^2,a = 

^^^M^ /^a^(pi)T„i(i^,2)/^a^(P2) i?:^^^"-^(^l) 
92 

+ X! / / dqWi^ hc,;^{p2)T„'^{Eg,Jha,'^{'K3{q'2q[x)) 

Go (tts (9291 a;), 92) Ha'^a,^{q'2q[x) {Tr4{q2q[x)q[a[\Ui^a 

- ^ J dx J dq2 q2'^ ^a;,(P2)Ta^(£'q^)^a;,(7r5(92'72 2;)) 
"2" 

Go(7r5(9292 a;), 92) I a'^a'^' i'l2q2 ^) (7^6(9292 a;)92 "2 1^^2,0 
= /^a'(P2)K.'(92) , (5.4) 



with 



K.'(9^) 



^^^2^r„,(i?,J/.„,(p2)i?:^^-^(0i) 

^2 



+ X! J J dq'iQi Ta'^{Eq'Jha'^{'K3{q'2q[x)) 

Go{T^3{qWlX),q'2) Ha'^a'.iqWl^) {T^4{q'2qiX)'lWl\Ul,a 
- ^ J dx J dq2 q2'^ Ta'^{Eq'^)ha'^{-K^{q'2q2 x)) 

"2" 

^0(71-5(9292 x)^l2) 4^03' (9292 x) (7r6(9292 x)q2 a-i \U2.a ■ (5.5) 
Then we insert the functions Z and Y under the integrals 

^o;(9i)- 

'^^'^^;'^^^ .',(9i)V,(pi)C";+"-(gi) 

+ y y d9k;'^T„;(i?,;)v,(^i(9W2a^)) 

Go(7ri(9i92a;):9i) G^^a^ ((?i92a;)/iai (7r2(9ig2a;))V'a^ (92): (5-6) 



and 



^^'^^;''^ a^(92)fe.>(P2)g:;,"'"-(gi) 

92 



^ J dx J dq[q^ Ta'^{Eq'^)ha[^_{TT'i{q2q[x)) 
Go (tTs (92 9^2;), 92) Ha,^a[iqWlX)ha'jTr4{q2qiX))Za[{q[) 

^ J dx J dq2 92 ^ Tai(-Bq^)/jQi (71-5 (9292 2;)) 
"2" 

Go (71-5 (92 92" > 92) ^a'a'" W2'l'2 X)h^''' Mq2<l2 x))V^''' {q'2 ) ■ (5.7) 



These two equations, Eqs. (j5.6l) and (j5.7p . form as set of coupled one-dimensional integral equations. The low order 
iterations exhibits the same features as discussed in detail in the previous section. Thus, we write 

Z^,{q[) = Z(?)(9l) + zW(9i) + + + • • • (5-8) 

and 

V^'M^) = V^Jil2) + V^l\q'2) + V^f{q'2) + V^f{q'2) + ■■■ . (5.9) 
From Eqs. (|5.6p and (|5.7p we can read off the different orders. For the lowest order we obtain, 

Z!^!i<l'i) = ^^^^^ (91) V,(pi) G;'f +"^(0i) (5.10) 
9i 



and 



V!^i<l^ = ^^^2_Mr„,(92)/..,(P2)i?:r-(^r) 

92 

- ^^^^^ r„^(92) ft„^(p2 )^:^"^(0i) . (5.11) 

92 



The second order is given by 
4?(9l) = 

(l + (-l)''i+'^'i) y y d9i92'EM(^9;) V,(^i(9i92a^)) Go(7ri(9'i92a^),9i) 
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and 



= (l + (-l/i+<) j dx j dq'^q'^^Y.^^'A^'iO Go{MQW2x),q[) 



'^Ga[a'^{QW2X) ha'^{n2{qW2X)) 

<5(«2 - 92) 



o r„.(,2)/..i(p2)i?^-"^(^i) 

'/2 



Q2 



= + ^^,(^;) y dx h^,{7ri{q[q2x))Go{MQl<l2x),q[) 

- Vi(7i"i(g'ig2a;)) Go{-Ki{(^^q2x),q[) 

Y^Ga>^a,_^{q[q2x) K'^{'K2{q'iq2x)) r^^fe) /la^fe) I?"^''"' (6*1) 



= (l + (-iyi+<) T^'Ml)T. 

y rfa; Vi(7i"i(g'ig'2a;)) Go{-Ki{q[q2x),q[) Ga'^a'^{q[q2x) K'^{-K2{q[q2x)) 

- T^'{h)K'{P2)Kr\e,) 



dx ha'^{ni{q[q2x)) Go{TTi{q[q2x),q[) Ga[a'2iq[q2x) ha'^{TT2iq[q2x)) 



J dx J dq'iqi^ Ta',^iq2) ha'^iTT3{q2q[x)) Go{-K3{q'2q[x),q2) 

- J dx j ^92'?^^ "^a^ (92) ^0^(7^5(929^2;)) ^0(775(9292' a;), 92) 

Z!^"ia;,"(9292"a;) (7^6 (9292" a;)) ^i°''(92") 

"2" 

= j dx j dq[q[^ Ta'^{q'2) K'^{n3{q2q[x))Go{TT3{q2q[x),q2) 

(9291a;) V, (7r4(gki^)) ^^^M ^^(p,) C7+"^(0i) 

- j dx j rfgag^^ '^0^(92) ^^0^(7^5(92922;)) Go(7r5(g292'a;),92) 



13 -^a^a;" (9292" a;) /i^" (^6 (9292" a^)) 

/// 
«2 

'^(92" - 92) 



^{q2 -92) 



9i 



Ta2"'(92) h^'''{p2) Dl 



-2 r„2"'(92)/i„;"(p2)I?:^?/'"^(^i) 

92 ^ 2 



y c^a; Tc^(g2) /la^ (71-3(92912;)) Go(7r3(g29ia;), 92) 
^ifa^ai(929ia;) V(7r4(929ia;)) ^(gi) V(pi) C";+"^(^i) 



(5.12) 
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III ^ 
Ta'^iq'l) ^a^(7I"5 (92922;)) Go{lT5{q'2q2X) , q'2) 

4" (9292.x) h^,,,ine{q'2q2x)) r„,"'(92) '"(^2) ^",?;'"^(^i) 
M(9^)Em(9i) C^;+"^^(^i) 

y rfx Vj, (71-3(92912^)) Go(7r3(g29ia;),92) H a.'.^a.'^'yW-vx) (7^4(922:)) 



.(«2)E 



Ta."'(92)/i '"(P2)i>::;?/'"^(^i) 



dx ha!^{Trs{q'2q2x)) 00(715(52922;), 92) ^a'^a'^' il'^l^x) (7r6 (92922;)) 

- r„.'"(92)/i "'(P2)i>::!?/'"^'(0i) 

y" dx /la^ (775(53922;)) 00(7:5(53922;), 92) 4^4" (92922;) 

V" (^6(92922;)) . 



(5.13) 



As for general forces the x-integration leads to logarithmic singularities in the external momenta. The next order, 
however, gives smooth functions 

^i?(^0 = 

+ r„i(50 JdxJ dq'2q'2^ KiiMqW2x)) Go{MqW2x),q[) 

J2Gaia'Mll2X) K'S-K2{q[(/2X)) V^^ {c/2) 

= (l + (-l)'i+«i) r^,{q[) j dq'2qi 

j dx K'^{ni{q[q2x)) Go{ni{q[q'2x),q[)'^Ga'^a!,{qW2x) K'^{n2{qW2x)) t„^(92) 

Ev/(9i) V(Pi)t^7+™^(^i) 

j dy /ia^ (773(53912/)) 00(773(53512/), 52) -ffa^ai (52912/) (774(922/)) 



-E 



Ta2"'(92) h^,,,{p2) oyr^e,) 



j dy /la^ (775(53522/)) 00(775(53532/), 52) -fc^c,;" (92922/) h^in{lTQ{^2<i2y)) 

ra,"'{Q2)h^iii{P2) Dyr^e,) 

(^2/ /lai (775 (92922/)) 00(775(53522/), 52) -^ala"' (92922/) ^ai" ('^S (9292^) ) 



(5.14) 



and 



^i?(92) 



y" d5i5iV„^ (52) /ic.^(773(525jx)) 00(773(535^0;), 52) 



12 



^ Ha'^a[{qWix) (71-4(92912;)) z^a'^ii'i) 

- y rfa; y" dq2 93 ^ 7-0^(92) K'^{T^b{q2q2 x)) 6*0(715(^292 a;), 92) 
Z]-^c<ic<;"(9292"a;) /i4"(7r6(9292"a;)) (92") 

/// 

= ■7'a;,(92) j dq'iq'i^ 

j dx /la^ (773 (929ia;)) Go(7r3(929la;),92)^i?aiai(929i2;) (71-4(92912;)) 



dy ha'^{TTi{q[q2y)) GQ{TTi{q[q2y),q[) G ^,^^n{q[q2y) h^,,{'K2{q'iq2y)) 
j dy /ia;(7ri(9'i92y)) Go(7ri(9;92y), 9i) G^,^^n{q[q2y) h^n{-K2{q'iq2y))] 

Toc'Ssll) j (^92" 92"^ 

/ dx ha'^{n5{q2q'2 x)) Go (775 (92 92" a:), 92) XI -^"^4" i'l^q-i x) {-K^ifi^q^ x)) 

HI 

(92") E '^"i (pi)^?™:^'"^ (^1) 

/i^" (7:3(92" 9iy)) GQ{'K3{q2qiy),q2) H ^'>> ^,^{({2 qiy) V (7:4(92"^)) 

- ^a;"(92')E r.,<<<<(92)/^„;'"(p2)^:^?r^(^^l) 

//// ^ 

j f^a'^' i'^5{q2 q2y)) Go (775 (92" 92^), 92") Ia'i'a'i"i^2 q2y) h^-" {■K6{q2 q2y)) 

- T„,""(92) V'(P2)£''^?"'"'(^l) 

2 a2 

dy hn, (7:5(92" 922/)) Go (775 (92" 92y), 92") -^a;"a"" (92"92y) 



/i'"' (7:5(92 92^)) 



(5.15) 



As shown in the previous section, there appear three-fold integrals. Two of them are over angles, leading to logarithmic 
singularities, and which arc then chminated by the third intermediate momentum integral. 
Thus we end up starting with n = 3: 

^S^('/i) - 

+ j dx j (i9292'M(9i) V,(7ri(9i92a;)) Go(7ri(9W2a;), 9i) 

EG„i„. (9'i9^a;) K'Mi<l'^^2x)) V^-'\<^2) 



Correspondingly, 



(l + (-l)''^+^9 r^'Sq'i) j dq2q2 J dx h^,{miqW2x)) Go{MqW2x),q[) 
Y,G^,^,{qW2x) h^,{MqW2^)) V^r'\q'2) • 

^i"^(92) = 



(5.16) 
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Y^H^'^^-iqWlX) K'S^,{q'^q'lx)) Z^^-'\q[) 



dx J dq^ ^ Ta'^{q2) ha'^i-K^iq'^q^ x)) Go(7r5(q^g2 2;), 92) 
y^^a'a':' (9292" h^''' {neiq'^q^ x)) ^'""^^(52") 



= Ta'^iq'2) J dq[qi J dx ha-jTT^iq^q'ix)) Go{7r3{q!2q'ix),q'2) 
^Ho,'^o^>^{q'2q[x) K'^{'K^{q'^q^x)) Z^^]^^\q^) 

- 7-0^(92) j dq2 q2^ j dx ha'^iiT^iq^q^ x)) G'o(7r5(g292 x),q'2) 

E 4" (9292" 2;) h.. {ii,{q'2q2 x)) vl;^,7^\q^) . (5.17) 

2 2 2 2 

/// 

Again, the singular integrals can be rewritten according to the method given in Appendix |B1 

VI. THE UNITARITY RELATIONS 

The scattering states depend on the initial state quantum numbers 

^(+) = ^(0) /g 

— Piqi,m2m3 V"--^/ 

like the initial state 

<I'0,a = <i>^,q„™,„3. (6.2) 

Using the full Green's operator, G = [E + ie — H)^^, to the full Hamiltonian, ^'p''iqi,m2m3 obeys the equation 

^ $" +GV^" (6 3) 

A second scattering state is defined by 

^(-) =$0 +G*V^° (6 4) 

*piqi,m2™3 Piqi,m2™3 ' Piqi,'n2™3 V"-^/ 

Both are related to each other as 

^(-) = +(G*-G)V^° 

*piqi,m2rn3 =^piqi,m2m3 ' \^ ^ I ' ^piqi, 77121713 

= *[,tqi,™2™3 + 2^*^(^ - ^)V^*p,q,,™2™3 ' (6-5) 

The S-matrix is defined as 

5''">3.™2m2 ^ ,^(-) |^(+) ^ .gg^ 

Piqi.Piqi — \ piqi, m^m^ I Piqi,ni2in3/- V"-"; 



Inserting Eq. (16. 5p leads to 



^117^171^,17121712 _ /!]/( + ) !*( + ) \ 

piq'i'Piii ^ p'iq'i>™2™3 piqi:™2"73' 
-27rz5(i?' - i?)($°iq;,™.,„. |V^|*(,t^^„,„3). (6.7) 



Now, we have due to general considerations 

>(t\ , ,14 , _ 

^ Pjq'j.mjmg I Piqi,m2"73/ \ p^q^.nijTiig I ^ piqi,m2m3 



(^d'i a', ,m'„m'„ 

(6.8) 



Consequently, 



C,77l2 1173, 7712 7112 _ I I^Q \ 

p'jq'j,piqi \ p[q[ ,m'2m'^\ piqi, 11721113' 
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-27ri5(i?'-£;)(<I>°,q,^„,„,|V|*W^,, 



~2mS{E' - i?)($;^;q-,„,.™. |[/"°Kq,.™,,„3). (6.9) 

For the last equation we used the definition of the transition operator J7°". 

Since the scattering states in the definition of 5" belong to the same Hamiltonian one has to have 

Viq'i,Piqi =Viqi,Piqi ^^i^-^): (O.iUj 

2 2 

where £' = ^ + (setting the a-particle as spectator), and correspondingly as similar expression for E' . The 
unitarity relation simply follows from the completeness relation spanned by the scattering states: 

\ pjq'j.mjm!, I piqi .mgnis/ 

7712 "^3 

or in terms of the 5*- matrix elements 

\ p^q j ,7712 Tjig I piqi ,m2 7?73/ 

E/ j3„" j3„" c™2"7n3",m2™3* C™2"ni3",ni2™3 ("fi 1 



'2 '"3 

rOO 



This can be rewritten in terms of the matrix elements of U . Using the completeness relation, 

y /dVi'dVi*" ) i = i> (6-13) 

Z ^ / ^ ^ ' p^ q^ ,7772 ^3 Pi qi '"^2 "^3 



and 



leads to 



8(E -i;')($°„ „ " "|(^°°l<'q' 777'7„') = 

($" „ „ „ „|<5(ffo-i^')f^"°l*p'a' 7„'777'), (6-14) 



+ 27r ^ J d^p'l(fq'i6{E" ~ E) 

7712 '^"'3 

^ Piqi,"J2™3 l^p'iqi,77l^7J7;j/ \^p"q",77l2 77I3 I I ^ p 1 qi ,77l2 TMa / 

= . (6.15) 
More interesting is the partial wave decomposed version for the on shell matrix element 

{pWia[\U°°\pmai) = U'J'^^{p[pi)5j,j,5M^M, (6.16) 



The dependence on p[,pi is sufficient since on shell qi — y 2Mi{E — ^),q'i — y 2Mi{E' ~ 
As shown in the Appendix [Cl the matrix element U^^'^^ (p'lPi) obeys the unitarity relation 

• rrOO,J*/ / \ ■ ttOO.J I I \ 

(Pl.Pl)-*C^aiai(Pl'Pl) 

+ 27rJ2j 5{E" i?)f/°?;J(p;\p;)C/°r (p;\pi) - 0. (6.17) 

The corresponding relation for the partial wave projected S-matrix element is 

y f SiE" - E)5{E" - E') ip'ip',) Si Jp[p,) = (6.18) 

„ J 1 ' 1 1 ' Pi Hi 

Note that not only discrete quantum numbers span the columns and rows of the S'-matrix but also the continuous 
quantum numbers p'lPi which describe how the energy is continuously distributed among the two relative motions. 
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VII. THE CAPTURE PROCESS N+N+o ^ '^HE 

The matrix element for the capture process is^ simply related to the time reversed photo disintegration process of 
^He into three free particles. It is well known [40] how to treat photodisintegration of ^He in the Faddeev scheme. In 
essentially the same manner one can formulate photodisintegration of ^He based on an effective three-particle picture. 
Let O be the photon absorption operator and jvl/eHe) the ^He ground state. The break up amplitude into nna can 
then be written as an infinite series of processes 

($0,a|[/o|*6Ho) = 

($0,a|O|«'«Hc) + 5]($0,a|K:GoO|*6He) 

i 

+ $]($0,a|V^^GoF,GoO|*6He>+--- . (7.1) 

Here Vi are the pair forces among the nn and na particles, and Go is the free propagator. This infinite series in terms 
of pair forces represents final state interactions (FSI). The first term is the direct break up process generated by O. 
Let us define 

($0,a|f/o|*6He) = (*0,a | O | ^-^He) + (*0,a | C/qz | *6He) , (7.2) 

i 

where C/gi comprises all terms with Vi to the very left: 

t/o^|*6He) = V,GoO\^m.) +V^Y, GoV, Go 0|*6He) + • • • . (7.3) 

j 

Clearly this can be summed up as 

C/0^|*6Hc) = V, Go 0|*6He) Go^C/o,|*6He). (7.4) 

3 

Separating the terms [7oi|^<5jj(,) to the left and introducing the i- matrices ti leads to three coupled Faddeev equations 
(1 = 1,2,3), 

[/o,|*6Hc) =t,GoO|*6He) +t»Go^[/o,|«'6Hc)- (7.5) 

The photon absorption operator O has to be symmetric under the exchange of the two neutrons, which wc number 
as particles 2 and 3. Thus, using the antisymmetry of |\E'6He) with respect to the two neutrons, one finds 

-P23C/02|*6Hc) = -f/03|*6He)- (7.6) 

This leads to the two coupled equations 

t/0l|^'«Hc) = tlGoO|*6He)+ilGo(l-P23f/o2)|*6Hc) 

[/02|«'6Hc) = i2GoO|«'6He) + i2Go(f/oi-P23f/o2)|*6He), (7.7) 

corresponding to Eqs. (|2.12p from Sect. [Ill 
The complete break up amplitude is given by 

($0,a|C/o|*6Hc) = (*0.a|O|^'6Hc) + ($0,a|f/0l|«'6Ho) 

+ (*0.a|(l-P23)f/02|*6He). (7.8) 

Using adequate pair forces and photon absorption operators (single particle currents, two-body currents and possibly 
beyond) these coupled equations can be solved by standard techniques (40j. 



VIII. SUMMARY 



The structure inherent in the continuum states of the n-|-n-|-a system has so far only been explored in the framework 
of the HH approach 0, H, [26l - l28i [3^ . There are strong initial and final state interaction peaks, not only in the nn 
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subsystem but also in the n — a subsystem. This poses a stih unsolved challenge for the expansion into the discrete 
set of K- harmonics as already known for the n + d^n + n + p system. This is pointed out by the authors of Ref. Q , 
who note, that even for a maximum K^ax — 20 in their calculation, the result is not completely converged. 

A corresponding investigation in the Faddeev approach is still missing. The aim of this paper is to lay the formal 
ground to do so. 

In the Faddeev approach all the structures in the relative motions of the three particles are mapped out correctly, 
thus leading to a reliable path to the three-to-three scattering S-matrix, which contains the information of the 
resonance structure of the ^He system. 

We derived two coupled Faddeev equations for the three-to-three scattering amplitudes. In a partial wave decom- 
posed representation they form a system of two-dimensional coupled equations for each fixed total angular momentum. 
The multiple scattering series being arranged in powers in the two-body t-matrices is the natural starting point for 
the solution of this coupled system of integral equations. The term linear in the t-matrices is disconnected. The next 
term, second order in t, has well established logarithmic singularities in the external momenta. Only the term of 
third order in i is a smooth function of the external momenta and thus can serve as driving term for the consecutive 
application of the Faddeev kernels. This provides all higher order terms which can then summed up by Pade. 

Since up to now nearly all Faddeev based investigations of the discrete structure of the n-|-n-|-a system are based 
on finite rank forces we also derived the continuum equations using this type of forces. The unitarity relations are 
especially interesting since the rows and columns of the S- matrix arc not only numbered by discrete quantum numbers 
but also by continuously varying on-shell momenta. 

Finally we provided Faddeev equations for the n+n+a capture process to the ^He ground state. We pointed out 
that it is not necessary here to first evaluate the three-to-three wave function as has been done in Ref. [llj, and 
that one directly can use a Faddeev form for the entire break-up amplitude with no disconnected terms as is done in 
modern calculations [s^, , and as it was pioneered in [sgI . [37| . 

This capture process is relevant for the production rate of ^He in astrophysical environments [4l| characterized by 
high neutron and alpha densities, e.g. those related to supernova shock fronts. In Ref. [2l[ this three-body process 
is approximated by sequential two-body processes, whereas in principle a genuine three-body reactions needs to be 
calculated. Very recently it has been pointed out [411 that currently employed two-step mechanisms over intermediate 
resonances in the three-to-three scattering of the n -\- n -\- a system are most likely insufficient, since the time delays 
for those intermediate steps are comparable to the duration of the entire process. This strongly supports the need for 
the approach we present in this paper. 
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Appendix A: Partial Wave Decomposition in the nna System 

Partial wave decomposition of three-body wave functions have often been documented, see for instance [39| and 
focus on the nna- system [46| . Thus we are here relatively brief. 

The projection of the free state, Eq. (|2.13l) onto the partial wave basis states of Eq. (|3.1I) is given by 



(pWiail^-piqimsma) = 

[l + ^(j^Ai J',^, M' - n){l[s'^j[, ^1 - 7712 - m3,m2 + TO3) 



11 , 

--51,77127713 



^{p'l -Pi) 5{q'i - qi) 



qf 



Y, 



(Al) 



This leads to 



(Pl^iail^ll^piqimams) = 



(^l^lii - m2 - 7773, 7772 + 777.3) 



1 1 

--S;^, 77127773 



Pi 



ql 



ql 



t:,{p,p,,E,,)Cy+"^^{p^q^), 



(A2) 



and gives the driving term of Eq. (13. 2p together with the amplitude C given in Eq. (13. 3p . 

In case of the second driving term in Eq. (j2.14p it is adequate to rewrite the free state |'l'piqim2m3) 

Jacobi momenta of the type 2 (where the neutron is the spectator): 



in terms of the 



P2 = -I3pi - 7qi 

q2 = Pi - aqi , 



(A3) 



with 



" = 2 



7 



777 + 771q 
2777 + 771a 

2(777 + 777a) 



(A4) 



and a/3 + 7=1. 
Then, 



I* 



Piqim2m3/ — 
|p2q2)|0)i|7772)2|m3)3 

lo(ii)i) , 



|p2q2)|0)i 17772)3 1^3)2 



(A5) 



where 



P2 = /3pi - 7qi 

q2 = -Pi - aqi . 



(A6) 



The partial wave projected state in system "2" is given by 



(P2'72a2l^piqim2m3) = 

S{P2 -P2) 5{q2 - 92) 



P2 



ql 
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^(P2 -P2) (5(92 - 92) . 

0..' 1 



~9 ~? 

P2 92 



^ij2l2J' , fJ-M' - A*)(^2^i2>A*- "I2,™2)yi^^_,„^(p2 



(A2^/2>^'^'-M-"^3,TO3)y;',A^,_^,_„^(q2) , (A7) 



and the second driving term becomes 

{pW2a'2\h\'^p,q, 

- ^-^^^^t^'Mp^,Eg,)Dy^"'%ep.,e^,), (AS) 

92 ^ 

with D and D given in Eqs. p.4p and p.Sp . For the kernel pieces we refer to (46| . 

Appendix B: Avoiding logarithmic singularities in the integrals 

We illustrate the new manner to rewrite the Faddeev kernel such that only a single pole singularity appears in an 
example (for more details see [H, E^) Consider the first kernel in Eq. p. 21) . of which the first piece can be rewritten 
as 

Jd.J dq'.qiY. J dp'ip? ^^^'^''"!2^'^'^"^^ t^MplE,0 Go{plq[) 
n ^ ' ' ^ /"^ ' '2 ^iP2-'^2iq'iq2x)) , , , 

Ga[a'^{qiq2x) / dp2P2 ^ (P2'72a2 I C^2,a ■ (Bl) 

J P2 

The two (5-functions are then changed according to 

// 

5{pi - '^i{q'iq2x)) = — rr 5{x - xo) Q{i - \xq\) 



S{P2 - T^2iqW2x)) = (5 1 P2 - \l iq? + 



a a 



with 



e ( 79? + ^P? - ^q2' ] , (B2) 
a a 



^0 = ^ — — ■ (B3) 



Inserting this into Eq. (jBTj) leads to 

J dx J dq'2qi^ j dpiPi^ ^^^^ (5(a; - xq) 6(1 - |a;o|) ta; (p'lPi , -E^j) 

Go(Pi,g'i) Ga'^a'^{q[q'2x) j dp'^p^ 6 - ^J"fq[^ + ^p^ - ^qi^ 
e ( ^q^ + ^-p^ - ^g;^) (p'2g^a^|C/2,. 
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Ga[a'AqW2^o){\hq? + -P? '72"2 



U2.a 



e 1 



"9 9 '9 '9 

Pi - a qi ~ 92 



a a 

© 1791^ + -Pi' ~ —qi 



2 , P "2 Pi '2 



(B4) 



The two 0-functions define the domain D for the integrations over and q'2 . Thus we end up with 

1 



[ dpip'i ta,{p'iPi,Eq,) 
aq'i J 



m 2 Ml 



/ c^'?292 V Ga'.a' ('j'l^ia^o) ( \/ 79i^ + -Pi q? q'20^2 



U2,a 



(B5) 



The singularity in Go(pi,<z'i) is now a single pole in pi for a given gj^. This type of singularity does not pose any 
numerical problem and can be implemented with standard techniques [50]. Note that for q'^ > \/1M\E there is no 
pole and one might as well keep the original form. 

Appendix C: Partial Wave Decomposed Transition Amplitude 

The definition of the partial wave decomposed transition amplitude is 

(*^'piqini2™3 1^ I$piqi)ri2m3) = 

^ y" dpiPl^rfgigi^ (^'p^q'^m^.m^ \pi il ^l) 

q"i 

(p'i'gl'c?i'|J7™|pigi(ii)(pigidi|$piqim2m3) . (CI) 

This inserted into Eq. (|6.15p yields 

8[E-E')^ ^Pl^l"l|*Pi1i'™2m3)(*p'iqi,mimyp-i'<fi'ai') 

I (pigiai|C/0"|pi'gl'aiO* -i(pi'gi'ai'|C/™b~igi«i) 
+ 27r ^ j d^p[d^q[ 5[e" - E) 

E / ('&1pi''9l''«i'V(Pi'''?i''ai''|C/™|Pi'9i'air 

Qi" 

J {^'' \pi''' qi" di'''){pi''' qi" di" \U°"\piqiai) 



Ql 

= 0, 



which, using the completeness relation, 

J2 J d^^p^d^qlSiE' -E){<S?''\p\''q\''di')A^''\p\'''q~i''di'') 

S{pi -Pi )5{qi -qi ) 



{p~i"r {qi"r 



leads to 



5{E - E')^! ^ j (Pigicfl|«'piqi,m2m3)($p;q;^m^m^|pi'gi'ai') 



(C2) 



(C3) 
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t{pmdi\U°"\pi'qi'aiT - *(pi'gi'ai'|C/"°|pigidi) 
2nY, I 5{E" -E){p\q^'d^'\U''%i'qi'd^r 



(pi qi di |[/°°|piqidi) 
. 



(C4) 



Then using Eqs. (|A1[) and (|6.16l) . the orthogonahty of the spherical harmonics and of the Clebsch-Gordan coefficients, 
one can project onto the on s-shell unitarity relation of Eq. ()6.17p . 
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FIG. 1: Diagrammatic representation of the first few terms of the multiple scattering series for the neutron-neutron-a system. 
Here the alpha particle (1) is indicated by the thicker line. 



